Scalings of domain wall energies in two dimensional Ising spin glasses 
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We study domain wall energies of two dimensional spin glasses. The scaling of these energies 
depends on the model's distribution of quenched random couplings, falling into three different classes. 
The first class is associated with the exponent 9 « —0.28, the other two classes have 9 — 0, as can 
be justified theoretically. In contrast to previous claims we find that 8 = does not indicate d — df 
but rather d < df , where di is the lower critical dimension. 

PACS numbers: 75.10.Nr, 75.40.Mg, 02.60.Pn 
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Spin glasses 0] exhibit many subtle phenomena such as 
diverging non-linear susceptibilities, aging and memory, 
making it a real challenge to understand these materials. 
In spite of much work, there is still no consensus even 
on the nature of the frozen order in equilibrium. More 
surprising still, the case of two dimensions also is not 
completely understood. In particular, the scaling of the 
stiffness, a cornerstone of spin glass theory, is different 
when the spin-spin couplings are of the form Jy = ±1 
compared to when they have a Gaussian distribution . 
This has been confirmed since using more powerful nu- 
merical techniques 0, 0] , and in fact it was interpreted 
in Q as a lack of universality, but this is unexpected and 
unexplained. Here we solve this puzzle: we find that dif- 
ferent types of quenched disorder lead to three distinct 
behaviors. In particular, we motivate why the class of 
models that includes the case Jy = ±1 gives for the stiff- 
ness exponent 9 = 0, and we explain what tells us about 
the lower critical dimension. 

The model, its properties and our methods - The 
model consists oi N = L 2 Ising spins Si = ±1 on a 
simple square lattice with periodic boundary conditions. 
The Hamiltonian is 



going from periodic to anti-periodic boundary conditions. 
At T = this reduces to measuring the difference 



H = — ^ Si J., 

(ij) 



ij Sj , 



(1) 



where the sum runs over all pairs of nearest neighbors 
(ij) and the Jij are the quenched random spin-spin cou- 
plings. We shall consider different distributions of these 
couplings, all of which are symmetrical about J = 0. We 
begin with continuous distributions; most common is the 
one where the Jij are Gaussian random variables with 
zero mean and unit variance. After that we investigate 
discrete distributions; the most common distribution of 
this type has Jy = ±1 with equal probability. 

An important feature of spin glass ordering is the spin 
glass stiffness; the corresponding exponent 9 describes 
how excitation free energies scale with the associated 
length scale. The standard way to measure this expo- 
nent is via the change in the system's free energy when 



6E = EP-E^ 



(2) 



where Eq P ^ and E^ AP%> are the ground state energies for 
the system with respectively periodic and anti-periodic 
boundary conditions say in the x direction. We are in- 
terested in the probability distribution of SE when con- 
sidering an ensemble of Jy and in the scaling law of its 
standard deviation AE: 



AE 



(3) 



Measurements of in two dimensional spin glasses (see 
for instance Q) give as —0.28. However, for the 
Jij = ±1 distribution, Hartmann and Young Q recently 
showed that AE remains of O(l) for increasing L, im- 
plying that in this case 9 as 0. In dimension d above the 
lower critical dimension d'f we have 9 > and spin glass 
ordering is stable against thermal fluctuations. On the 
contrary, when 9 < 0, thermal fluctuations prevent spin 
glass ordering. Because of this, the authors of conjec- 
tured that di = 2 for the Jij — ±1 model. We shall see 
that d\ should be identified with the highest value of d 
where 9 < 0, and so in fact df as 2.5 as believed before 
the study in Q. 

In this work we address these questions by first deter- 
mining numerically the properties of P(SE) and then by 
using the real space renormalization group picture. For 
the first part, we compute the ground states of our sys- 
tems using a heuristic algorithm In practice, when 
the lattice is not too large (L < 80), the algorithm re- 
turns the ground state with a high level of confidence 
for all of the distributions we shall consider in this work. 
The problem is to reduce enough the statistical errors; in 
practice we used a few tens of thousands of samples at a 
few values of L for each case. 

Class 1: "continuous" distributions - We first fo- 
cus on distributions P(J) that include a continuous part 
(we shall see later that this class includes certain discrete 
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distributions also). When L is sufficiently large, SE can 
then take on arbitrary values. The value of 9 for contin- 
uous distributions is well known only for Gaussian J^ ; 
in fact we are aware of no tests of universality in d = 2, 
though the standard lore is that both 9 and the shape of 
P(SE) are universal Q|. 

In a first series of runs we obtained P(8E) and AE for 
the model with Gaussian couplings. Then we moved on 
to a continuous yet singular probability density P(Jy): 
P(Jij =J) = f/i(J) + (1 - /) W)> where P X (J) 

= [e^+e^J/VSi, P 2 (J) = [<y(J-l)+*(J+l)]/2, 
and / is a measure of the height of the distribution at J rs 
0. We refer to this P(Jij) as the broadened bimodal (BB) 
distribution since it reduces to the Jy = ±1 distribution 
when / = 0. 

In Fig. ^ we show AE function of L when 

P(Jij) is: (1) a Gaussian of zero mean and unit variance 
(GAUSS data); (2) the BB distribution, with / = 0.1 
(BB 0.1 data); (3) as in (2) but with / = 0.2 (BB 0.2 
data); (4) Gaussian but with the part in the interval 
[-0.5,0.5] forced to be (HOLE data). Note that this 
last distribution has a large gap around Jij = 0. In the 
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FIG. 1: AE as a function of the system size for four different 
Jij distributions. Straight lines are best one-parameter fits of 
the form const x L -0,282 . Inset: the probability distribution 
P(\5E\/AE) at L = 40 for three of these distributions. 

Gaussian case the power law scaling of AE can be de- 
termined with good accuracy already from quite small 
lattices; fits to these data lead to 9 = -0.282 ± 0.004, in 
agreement with previous work. The distributions (2), (3) 
and (4) give rise to a similar scaling albeit only at larger 
L values. We have also considered other distributions 
such as P(Jij) uniform in [—1.5, —0.5] U [0.5, 1.5] (notice 
that this distribution also has a gap around J = 0), ob- 
taining similar results. It thus seems very reasonable to 
expect that all distributions with a continuous part will 
lead to the same exponent, 9 as —0.28. 

A second universality issue concerns the shape of 
P(SE). In the inset of Fig. we show the probabil- 



ity density P(\SE\/AE) when L — 60 for the BB 0.1, 
BB 0.2, and GAUSS data: the different data sets basi- 
cally coincide within statistical errors, strengthening the 
claim that in this class the distribution of domain wall 
energies is universal (the curve displayed is just to guide 
the eye). 
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FIG. 2: AE/Q (top) and P(SE = 0) (bottom) as a function of 
the system size for three discrete Jij distributions: ±1 (PM1), 
diluted ±1 (DIL), and ±1, ±2 (PM1PM2). 

Class 2: quantized energies - At variance with the 
former distributions, the Jij = ±1 model leads to 9 ~ 
0. We show in Fig. EJthat in this model AE saturates 
quickly as L grows. Is the Jij = ±1 model a special case, 
a class on its own? The crucial point is that the possi- 
ble SE values are quantized: SE is always a multiple of a 
quantum Q, here Q = 4. This led us to consider distribu- 
tions other than the ±1 one with this same quantization 
property. We begin by "diluting" the = ±1 model, 
setting J^ = with probability 0.2. The main effect of 
this is to reduce the quantum from 4 to 2; indeed, the 
local fields now can take the value 0,1,2,3,4 instead of 
0,2,4. In Fig. 0we see that for this model (DIL data) 
AE seems to saturate, so again 9 = 0. However the con- 
vergence is slow. In any renormalization group picture 
this convergence is governed by a "correction to scaling" 
exponent lu. We assume 9 = and that the asymptotic 
value of AE is a non-zero constant given by the J = ±1 
data; then we fit the diluted (DIL) model to the form: 

AE(L) fa AE(L = oo) + AL~ U . (4) 

with A and u> adjustable parameters. We have also con- 
sidered distributions where Jij = ± J\ or ± J2 with equal 
probability (we have studied the cases J2/J1 = 1.5, 2 and 
3.). Again we find the convergence to be slow but fits as 
in work well; furthermore, all the estimates of lu are 
similar, being in the [0.4 : 0.6] interval. All these facts 
justify the claim that 9 = whenever SE is quantized. 

Just as in the continuous case, to analyze the shape of 
the distribution of 8E we must choose a scale; the cor- 
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rect choice is to compare the histograms after measuring 
all energies in units of the basic quantum Q. To test 
whether the histograms for the different Jy distributions 
become identical in the large L limit we plot in Fig. |3 
(lower panel) the probability P(5E — 0) to find a zero en- 
ergy domain wall. The data suggest that the histograms 
become identical in the large L limit, i.e., they support 
universality. (Following Eq. we fix the asymptotic 
value of P(SE = 0) to be that given by the J = ±1 
model, and then we determine u>; in the plot we show 
these fits; they are all good and the values of u are close 
to 0.5.) We have checked in detail that this claim applies 
to the quantized distributions mentioned before and to 
the DIL model with 10% dilution. 

Class 3: quantized energies revisited - So far we 
have only considered situations with even values of L. If 
L is odd (and the = ±1), the possible values of SE/ J 
are ±2, ±6, ±10, ...The quantum Q is still the separa- 
tion between the energy values, but the positions of the 
histogram entries are different (in particular, 5E = is 
not allowed) . A somewhat trivial consequence of this is 
that necessarily 8 > as AE/Q is greater or equal to 1/2 
for all L. Consider now the question of the universality 
of the histograms. We have checked within our error bars 
that the large L limit of P(SE/Q) for the Jy = ±1 model 
is the same as that obtained using the J2/J1 =3 model 
(still with L odd of course). This kind of quantization 
thus gives rise to a third class, again with 9 = 0. 

Could there be further classes with quantized energies? 
Since we have imposed reflection symmetry of the distri- 
bution of the the only possible histograms are the two 
we discussed: SE is a multiple of the quantum Q or of the 
form (n + l/2)Q, where n is integer. If the universality 
class depends only on the possible histogram types, then 
no other classes arise. 

Discrete does not mean quantized — Let us also con- 
sider the case where the couplings are discrete but where 
there is no quantization. We consider the distribution 
P(J) = \[5{J ± Ji) + 5(J ± J 2 )] (IRR for "irrational" 

hereafter), where J\ = 1 and J2 = 1+ 2 ~ 1-618 is the 
golden mean. Clearly we have SE — 1(nJ\ + mJ^) where 
n and m are integers. Since J2/J1 is irrational, the set 
of possible SE values becomes dense when L — ► 00 and 
so it is natural to conjecture that this P(Jij) leads to do- 
main wall energies in class 1. Our findings are that AE 
decreases with L and shows no sign of saturation, and a 
power law fit gives 8 = —0.29 ± 0.01, the value associ- 
ated with class 1. Our conjecture is thus substantiated 
by these findings. 

The convergence of P(SE / AE) to its limit is more 
problematic: for finite L, the distribution is the sum of 
a finite number of delta functions: we can only hope 
to have a "weak convergence" to the P obtained with 
the Gaussian couplings. In these conditions it is appro- 
priate to consider the integrated probability distribution 




SE/AE 

FIG. 3: The integrated probability distribution for SE/AE in 
the Gaussian case, in the two BB cases, and in the irrational 
J2/J1 case (L — 60). Inset: binned probability distribution 
of SE/AE for the irrational case (IRR), L = 60. 

n<(X) = Jip(^)i(^). In Fig. Ewe plot n< for 
the irrational and for some continuous cases. The plots 
are very similar, supporting the claim that the discrete 
distribution IRR leads to domain wall energies in class 1 . 
We have also included P{8E / AE) as an inset into Fig. 
we have used a small bin size that allows one to observe 
the complex structure. 

The case of hierarchical lattices — The effect of hav- 
ing quantized SE can also be studied on hierarchical lat- 
tices. One advantage is that one can study very large 
sizes, a second is that one can access a continuous range 
of dimensions. We have focused on Migdal-Kadanoff lat- 
tices [6j; these are obtained by recursively "expanding" 
graphs. Starting with one edge connecting two sites, one 
replaces it by b paths in parallel, each composed of s 
edges in series, leading to b ■ s new edges. This procedure 
is repeated hierarchically; after G "generations" the dis- 
tance between the outer-most spins is L — s G , while the 
number of edges of the lattice is (b ■ s) G . The dimension 
of these lattices is d = 1 + ln(6) / m(s) . One puts an Ising 
spin on each site and a coupling Jy on each edge. Peri- 
odic boundary conditions simply imply that the two end 
spins must have the same value, and from this we define 
SE. The probability distribution of SE can be followed 
from G to G + 1. The recursion equations for P(SE) 
make sense for any s pj: s can be an integer but it can 
also be any positive real value! One may then compute 8 
for an interval of dimensions, usingeither continuous Jy 
(for instance to check universality |a,|9() or quantized Jy 
couplings (our focus here). 

In Fig. 21 we show 8 as a function of dimension d (s is 
variable, b is fixed and set to 3). We show the values for 
continuous distributions and for when the quantization is 
of the form of class 3. As expected, if in one class 6 > 0, 
all classes lead to the same value of 8, i.e., quantization 
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is irrelevant when the energy scale diverges. However, as 
soon as 9 < in the continuous case, quantization gives 
rise to a histogram fixed point distribution in which the 
\SE\ are concentrated on the few lowest values and 9 = 0. 

Similar results are obtained for class 2 quantization 
but there is an interesting difference. Indeed, since SE 
can be zero in this class, one sees two further fixed points. 
An obvious one is associated with having P(5E = 0) = 1, 
i.e,. all domain wall energies vanish. It is easy to see that 
this fixed point is stable and has 9 = — oo; there is no 
spin glass stiffness, and the system is paramagnetic even 
at zero temperature. The other fixed point is unstable 
and has 9 = 0. What is the interpretation of these two 
extra fixed points? To allow SE to be zero, one can 
think of the diluted model where some of the bonds have 
Jij = 0. Clearly when the dilution is strong enough, the 
non-zero bonds will no longer percolate and we are in a 
strongly paramagnetic phase; the renormalization group 
(RG) flow in this phase takes one to the P(SE = 0) = 1 
fixed point. On the contrary, at low dilution, we are 
in a spin glass phase and the RG flows are towards the 
other stable fixed point. On the boundary of these two 
phases, the RG flows takes one to another fixed point 
which is unstable: it is associated with the paramagnetic 
to spin glass transition as dilution is decreased. Such 
considerations have previously been developed for d = 3 
Migdal-Kadanoff lattices . 

Finally, we see that it is appropriate to define the lower 
critical dimension df from the end point of the 9 = 
curve; 9 = on its own does not signal d = df. 
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FIG. 4: 9 as a function of d for a one-parameter family of 
Migdal-Kadanoff lattices; we display two sets of data points, 
one for continuous Jij distributions, the other for quantized 
distributions. 

Discussion — Our numerical evidence of universality 
(both for 9 and P{SE / AE)) is very strong for continuous 
and related distributions (see figures ^ an( l El • But we 
also find universality classes when SE is quantized. This 



classification is substantiated by the behavior of 9 and 
of the fixed point distributions of domain wall energies 
in Migdal-Kadanoff lattices. It is appropriate however to 
be cautious and to remark that the correction to scaling 
exponent u> we measure (see equation^} is small, u) ~ 0.5. 
Because of that we are not able to completely exclude 
the Bray and Moore expectation that AE(L = oo) = 
[lll | . Our most extensive data are for the model DIL with 
/ = 0.2. Here our fits give AE(L = oo)/Q = 0.49(1) 
while if we force AE(L = oo) = 0, the % 2 of the fit 
increases by 2.3; thus AE(L = oo) = is not excluded 
by our data though it appears as much less likely. 

What is the source of the universality we observe? In 
the Migdal-Kadanoff lattices, the renormalization group 
transformation is clear and so the different classes are 
very natural. For the Euclidean lattices the existence of 
a renormalization group transformation for SE has not 
been established, but since our data point to universal- 
ity, it should be possible to define such a transformation. 
Note that its fixed point (and thus P(SE/AE)) will de- 
pend on the aspect ratio and on the fact that we use pe- 
riodic boundary conditions. Our P(5E) are thus a priori 
not comparable to those of 4] where one of the directions 
had free boundary conditions. 
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